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ABSTRACT 

Threshold enhanced perturbative QCD corrections to rapidity distributions of Z and W ± bosons 
at hadron colliders are presented using the Sudakov resummed cross sections at N 3 LO level. We 
have used renormalisation group invariance and the mass factorisation theorem that these hard 
scattering cross sections satisfy to construct the QCD amplitudes. We show that these higher order 
threshold QCD corrections stabilise the theoretical predictions for vector boson production at the 
LHC under variations of both renormalisation and factorisation scales. 
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Recent theoretical advances in the computations of higher order radiative corrections in pertur- 
bative Quantum Chromodynamics (pQCD) have lead to extremely accurate predictions for several 
important observables needed for physics studies at the Tevatron collider in Fermilab as well as 
at the upcoming Large Hadron Collider (LHC) in CERN [1]. The Drell-Yan (DY) production of 
di-leptons [2], which is one of the dominant production processes at hadron colliders, can be used 
to precisely calibrate the experimental detectors. In addition, the DY process provides precise 
measurements of various standard model parameters through measurements of the rapidity dis- 
tributions of Z bosons [3] and charge asymmetries of leptons coming from W boson decays [4]. 
Possible excess events in di-lepton invariant mass distributions can point to physics beyond the 
standard model, such as R -parity violating supersymmetric models, models with Z', or with con- 
tact interactions [5] and gravity mediated models (see [6] for recent update). The precise measure- 
ments of Z and W boson production cross sections, various distributions and asymmetries by both 
the DO and CDF collaborations [7] at the Fermilab Tevatron, where y/S = 1 .96 TeV, have already 
provided stringent tests of the standard model. These have already played an important role in 
bounding the mass of the Higgs boson. Similar measurements at the LHC will provide even more 
stringent tests due to the increase in the number of events at a/S =14 TeV. 

The total cross sections for the Z and W ± production are known in pQCD up to next-to-next- 
to-leading order (N 2 LO) [8-14]. Resummation programs for the threshold corrections to the total 
cross sections for DY production of di-leptons are also known [15, 16] (see also [17]) and one 
can consult [18, 19] for next-to-next-to-leading logarithmic (N 2 LL) resummation results so it is 
straightforward to study the threshold effects in Z and W ± production at the cross section level. 
Recent QCD results at the three loop level [20-26] have lead to predictions for the resummation up 
to N 3 LL [27-30]. Notice that the fixed order partial-soft-plus-virtual N 3 LO corrections [27,30] to 
the Higgs and DY production show the reliability of the perturbation theory results and demonstrate 
stability against the variations of renormalization and mass factorization scales. Exact results up 
to N 2 LO are also available for less inclusive observables for di-lepton [31], Z and W ± [32] (see 
[33-35]) production. Recently the dominant QCD threshold corrections to the rapidity distribution 
of di-leptons in the DY process at N 3 LO have been obtained in [36]. It was found that these 
corrections are indeed small and reduce the scale uncertainties significantly making the predictions 
more reliable. The fixed order results as well as the resummed results reveal very interesting 
structures in the perturbative QCD series (see, [37-42]). 

The hard scattering cross sections computed using the QCD improved parton model are often 
sensitive to variations in the renormalisation and factorisation scales usually denoted by /jr and /jf 
respectively. The former originate from ultraviolet renormalisation while the latters originate in the 
mass factorisation of collinear singularities. In addition to the scale uncertainties the fixed order 
computations suffer from the presence of various large logarithms which arise in some kinematical 
regions. These regions are often important from the experimental point of view and these large 
lorgarithms, which spoil the standard perturbative predictions, should be resummed in a closed 
form. For instance resummation formulae supplemented with fixed order results can predict the 
dominant higher order threshold corrections to various observables. These threshold corrections 
are large when the fluxes of the incoming partons are large, which occurs at large partonic energies. 
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In [37] we computed the soft distribution functions that resum the soft gluons coming from 
real gluon emission processes in DY production and Higgs production and also found that they 
are related by the colour factor Ca/Cf, see also [43] and [44]. Using the soft distribution func- 
tions extracted from DY, and the form factor of the Yukawa coupling of Higgs to bottom quarks, 
we predicted the soft-plus-virtual (sv) parts of the Higgs production cross section through bottom 
quark annihilation beyond N 2 LO with the same accuracy that the DY process and the gluon fu- 
sion to Higgs process are known [27,30]. This approach was then successfully applied in [38] 
to Higgs decay to bottom quarks and hadroproduction in e + e annihilation. Since our results 
in [30, 37, 38] are related to that of the standard threshold resummation, we could determine [37] 
the threshold exponents D\ up to three loop level for DY and Higgs production using our re- 
summed soft distribution functions and the quantities B\ for deep inelastic scattering, Higgs decay 
and the hadroproduction of Higgs bosons. In [36] we extended this approach to include xp and 
rapidity differential cross sections for di-lepton pairs in DY production and for Higgs bosons in 
Higgs production processes. In this paper we apply these same methods to study the effects of the 
dominant threshold corrections at N 3 LO to the rapidity distributions of the Z and W ± bosons in 
hadron-hadron collisions. 

In [36] we formulated a framework to resum the dominant soft gluon contributions coming 
from the threshold region to the xp and rapidity distributions of DY di-lepton pairs and Higgs 
bosons at hadron colliders in the zi (i =1,2) space of the kinematic variables. We recapitulate the 
main points here to make this paper more understandable. The threshold region corresponds to 
Zi — > 1 and in this region all the partonic cross sections are symmetric in z\ <-> zi- To obtain the 
resummed result, we used renormalisation group (RG) invariance, mass factorisation and Sudakov 
resummation of QCD amplitudes. Using the resummed results in z% space we predicted the sv parts 
(also called threshold corrections) of the dominant partonic xp and rapidity distributions beyond 
N 2 LO. We follow the similar approach here to obtain the dominant threshold corrections at N 3 LO 
level for the rapidity distributions of Z and W ± bosons at both the LHC and Tevatron energies. 
See [45] for an early reference where the resummation for DY differential distributions at rapidity 
Y = (or xp — 0) was considered. 

The differential cross section for producing a vector boson can be expressed as: 

d2 ° J <jxUW)w\xUW), (i) 



dq 2 dy 

where q is the four-momentum of the vector boson. In our case q 2 = M 2 where J = Z, W ± but 
for convenience we use q 2 for most of this paper. Later we will present plots for d 2 o J /dqdy, 
where q now represents ^/q 2 = Mz or My/ for Z and W ± respectively. Our normalisation is 
W^ om {x\,x^_,q 2 ) = 8(1 — jq)8(l — JGj). The superscript / represents light-quarks (q), gluons (g) 
and heavy quarks (b) but we only need the first case here so / = q for the rest of the paper. The 
xf (i = 1,2) are related to q 2 , the scaling variable x = q 2 /S, and the rapidity y of the vector boson 
J: 

1 ( P2-q\ i f x°A T _ 00 
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Here S = (pi + P2) 2 is the square of the hadronic center of mass energy and p, are the momenta of 
incoming hadrons Pj (i = 1,2). 

In the QCD improved parton model, the function W^x^x®, q 2 ) can be expressed in terms of the 
fitted parton distribution functions (PDFs) appropriately convoluted with perturbatively calculable 
partonic differential cross sections denoted by A ! d ah as follows 



W l {x\,x\,q 2 ) = £ f dx x [ dx 2 ^ b {x h x 2 ,M 2 F 



X 



f 1 f l 

' dzi dz 2 5(xj-xizi) 8(^2-^222) Aj flfe (zi, Z2,q 2 ,MF^l), ( 3 ) 
Jo 



where the subscript d denotes the particular differential distribution one is studying (y, xf etc). 
Here hr is the renormalisation scale and hf the factorisation scale. The function ^^ b (xi,X2 : p F ) is 
the product of PDFs f a (xi : p F ) and fb(x2,M F ) renormalised at the factorisation scale /j f . That is, 

X« b (x h X 2 A) = fa l {xiA)fl 2 {x2A), (4) 

with Xi (i =1,2) the momentum fractions of the partons in the incoming hadrons. 

The partonic cross sections can be expressed in terms of soft and hard parts. The soft parts 
come from the soft gluons that appear in real emission as well as in the virtual processes. The 
infra-red safe contributions from the soft gluons can be obtained by adding the soft parts of the 
differential cross sections with the ultraviolet renormalised virtual contributions and performing 
mass factorisation using appropriate counter terms. These combinations are called the "soft-plus- 
virtual" (sv) parts of the differential cross sections. Hence we write 

^d,ab{zhZ2,q ,/J F ,(J R ) =A I a ab(zhZ2,q ,PFiHR)+Kb^dM^ Z2 ^ »i"F»i"j?)> l = c l- ( 5 ) 

The hard parts of the differential cross sections A^ d (zi,z2,^ 2 ,A , F'A'i) can be obtained by the 
standard procedure(see [6,46]). The sv parts of the differential cross sections are obtained using 
the method discussed in the [36] so that 



A s ^(zi,Z2,4 2 ,a4a<f) = Cexp (v I d (q 2 ,/j 2 R ,iJ Z F ,zi,Z2,z)^ 



(6) 

£=0 



where the x ¥ 1 d (q z ,fi R ,y: F7 zi,Z2,£) are finite distributions computed in 4 + £ dimensions and they 
take the form 

^(<7 2 ,AiA^i,Q,e) = (ln|F / (a,,e 2 , A / 2 ,£)| 2 ^5(l-zi)5(l-z 2 ) 

+2 <Z> d I (a s ,q 2 ,iJ 2 ,z h Z2,£) - C\nr n (a s ^ 2 ^ 2 F ,zu£) 5(1 -z 2 ) 
-ClnT n (a s ,iJ 2 ,iJ 2 F ,Z2,e) 5(1 -zi). (7) 



3 



The symbol "C" means convolution. For example, C acting on the exponential of a function 
f(z\,zi) means the following expansion: 

Ce f(zi,z 2 ) = 8 (i_ Zl ) 8 (i_ Z2 ) + i. / ( Zl>Z2 ) + i. / ( Zl)Z2 ) (g)/ ( Zl>Z2 ) 

1 



+ gy/(zi»Z2)®/(zi,Z2) 



)/(zi,z 2 ) + ---. (8) 
In the rest of the paper the function f(zi,Z2) is a distribution of the kind 8(1 —Zj) or Di(zj), where 

"hV'(l-z,- 



*>i(zj) = 



i = 0,l,- 



and ;'=1,2, 



(9) 



and the symbol ® means the "double" Mellin convolution with respect to the variables z\ and zi- 
We drop all the regular functions that result from these convolutions when defining the sv part 
of the cross sections. The F / (a i .,2 2 ,/j 2 ,e) are the standard form factors coming from the purely 
virtual parts of the cross sections. In the form factors, we have Q 2 = —Mj. The partonic cross 
sections depend on two scaling variables z\ and zi- The functions 4> d / (a J ,^ 2 ,/j 2 ,zi,Z2,e) are called 
the soft distribution functions. The unrenormalised (bare) strong coupling constant a s is defined as 



a x = 



~2 
Ss 

\6% 2 



(10) 



where g s is the strong coupling constant which is dimensionless in n = 4 + £ space time dimensions. 
The scale /j comes from dimensional regularisation which makes the bare coupling constant g s 
dimensionless in n dimensions. The bare coupling constant d s is related to renormalised one by 
the following relation: 



Sza s =Z{i? R )a s {}? R ) 



(ID 



where 5e = exp { | [y E — ln47t] } is the spherical factor characteristic of n-dimensional regularisa- 
tion. The renormalisation constant Z(//|) relates the bare coupling constant d s to the renormalised 
one a s {j^). They are both expressed in terms of the perturbatively calculable coefficients p\ which 
are known up to four-loop level [47,48] in terms of the colour factors of SU(N) gauge group: 



C A =N, 



C f = 



N 2 -l 
IN 



1 



(12) 



Also we use n f for the number of active flavours. 

In dimensional regularisation, the bare form factors F 1 (d Sl Q 2 ,/j 2 ,e) satisfy the following dif- 
ferential equation [49-52]. 



e 2 ^lnF'(a,,eV,e) 



a ^ £ 



+ G 1 



a 



Si 9 j 7 : 

4 v 2 



(13) 
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The fact that the F 1 (d s ,Q 2 ,/j 2 ,£.) are renormalisation group invariant and the functions G 1 are 
finite implies that the K 1 terms can be expressed in terms of finite constants A 1 , the so-called cusp 
anomalous dimensions and the coefficients p\. The formal solution to the eqn. ([T3l . in dimensional 
regularisation, up to four- loop level is obtained in [22,37]. The finite constants G/(e) (see eqn.(19) 
of [30]) are also known [23] to the required accuracy in e. These constants G/(e) are expressed 
in terms of the functions B\ and //. The B\ are known up to order a] through the three-loop 
anomalous dimensions (or splitting functions) [20, 21] and are found to be flavour independent, 
that is B q = Bj . The constants fj are analogous to the cusp anomalous dimensions A' that enter 
the form factors with A q = A\. It was first noticed in [43] that the single pole terms in £ in the 
logarithms of the quark and gluon form factors up to two-loop level (a 2 ) can be predicted by the 
Cf — > Ca substitution. The structure of single pole terms of four-point amplitudes at the two-loop 
level can be found in [53, 54]. The UV divergences present in the form factor are removed when 
the bare coupling constant d s undergoes renormalisation via the eqn.(fTTT). 

The collinear singularities that arise due to the presence of massless partons are removed using 
the mass factorisation kernels T(zj : /J F: e) in the MS scheme (see eqn.([7J)). We suppress their de- 
pendence on d s and jj 2 . The factorisation kernels r(zj,/u F ,£) satisfy the following renormalisation 
group equations: 

^ F ^T{zj^ 2 F ,z) = \p{zj^ F )®T{zj^ 2 F ,z) 7 = 1,2, (14) 

where the P(zj,fJ F ) are the DGLAP matrix-valued splitting functions which are known up to three- 
loop level [20,21]: 

oo 

p(z J A) = J L" i sU)p {i - l) (z J ). d5) 

i=\ 

The diagonal terms in the splitting functions P^ l \zj) have the following structure 



pinzj) 



B! +1 5(l-z j )+A! +l 'D (z j ) 



+ PrL(Zj), (16) 



where jj(zj) are regular when the argument approaches the kinematic limit (here Zj — > 1). The 
RG equations can be solved by expanding them in powers of the strong coupling constant. Only 
the diagonal parts of the kernels contribute to the sv parts of the differential cross sections. We find 
the solutions contain only poles in £ in the MS scheme: 

r(zjA,*)=W-zj) + i*,(^) T 4r (0 (zy,e). (17) 

An expansion for the rW (zj, e) in negative powers of £ up to four- loop level can be found in [37]. 
The rn(d s ,iLi 2 7 jJ F ,Zj,£) in eqn.© is the diagonal element of r(zj,/d F ,e). 

From the eqn.(fT3l) and the fact that the A^ v l are finite in the limit £^0we obtain 



K d [ d s ,-p-,ziiZ2,£J +G d ^a i .,-2,-2,zi,z2,e 



(18) 
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where now the constants K d contain all the singular terms in £ and the gJ are finite functions of 
£. The functions <t>j(d s ,q 2 ,jj 2 ,zi : Z2,£) also satisfy the renormalisation group equations: 

/4-^2<&M,$V,zi,z 2 ,e) = 0. (19) 

The <t> ( j(d s ,q 2 ,fj 2 ,z\ 1 Z2,£) should contain the correct poles to cancel the poles coming from F 7 ,Z 7 
and Tji in order to make A s ^ l finite. This requirement unambiguously determines all the poles 
of this distribution. The solution to the Sudakov differential equation for the soft distribution 
functions in eqn.(fT8l) can be written as 



i=l 

where 



(21) 



The constants (e) are expanded in powers of the bare coupling constant d s as follows 



.£ 

(2 \ oo / 2 \ 'T 

d s ,^,z h Z2,Ej =5(l-zi)8(l- Z2 )X;4^J 5^Zi' (i) (£). (22) 

Using the RG equation for K d (d Sy f^ ( /iJ 2 ,zi,Z2^), one finds that the constants K^'\e) are iden- 
tical to K I,{i \e) given in [30]. The constants g/' (e) are related to the finite boundary functions 
G rf (a s (^ 2 ),l,zi,z 2 ,£). We define the Q di {£) through the relation 

.£ 

I^ K "^ 1 "^ )' 7 ^ ^) = I^(^ 2 (l-zi)(l-z 2 ))^(e). (23) 
i=l V M / i=\ 



i=l V i" / /_i 
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We obtain the z\,Z2 independent constants Q d -(e) by demanding the finiteness of A s ^ given in 
eqn.©. Before setting e = in eqn.©, we expand A s ^ l as 



& S d,l(ZhZ2,q 2 ,l4,l4,z) = Y, ai s(/4) AS df ] {zuZ2,q 2 ^l^ 2 F ^) ■ (24) 

(=0 



Using the above expansion and eqn.© we determine these constants using the known information 
on the form factors, the mass factorisation kernels and the coefficient functions A^'j 1 ^ expanded 
in powers of £. The structure of the GJ(e) in the form factors involving the constants f and p\ was 
given in [30]. The constants Q in the soft distribution functions also have a similar structure: 



ill® = -f(+L^;i 
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i 

k=l 



k=\ 

OO 

?i(e) = -A ' 2Pi ?if" - 2Po (eii" + 2P„? if) + £ e^/f , 

k=\ 

9d,M = -/l-2P25i; 1 (1) -2|3 1 (5i i f +4ft,7/f'), 

-apo^f'+apo^'i 2 * +4^« 3 ' ) + £ e's/f . (25) 

The terms proportional to £ at every order in a s are determined using the known o J to order N 2 LO 
and the following identity: 

! l dx° x f l dx° 2 (x^f- 1 ^L=[ 1 dz x N - { a J . (26) 
Jo Jo ay Jo 

We find 

F <7,(2) _ r (K\ 

J d q f = c F (^g) , 

( 328 10„ 8„ \ 
+C F n f ^- — + ^ 2 + -^. (27) 

Using the resummed result given in eqn.©, the exponents g ; 9 (£) (see [23]) and Q / ? -(e) , we can 
obtain the higher order sv contributions to the differential cross sections. The available exponents 
are 





7r?.(y) 


for 


; = all 




7T 9)0') 

y d,i 


for 


7 = 0,1 




?d,3 


for 


7 = 



in addition to the known (3/ (z = 0, 1,2, 3) and the constants in the splitting functions A?, 5f (i = 
1,2,3) and jf (i = 1,2,3). The constants g% J are known for j = 2,3 also (see [22]). Using our 
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approach we have obtained the exact A^ v ^ up to N 2 LO (z = 0, 1,2) [32]. The coefficient of the 

8(1 — zi)8(l — zi) part depends on the constants g\ jgf' 1 , q\ which are still unknown for 
N 3 LO, so we can only obtain a partial result for A , , i.e., aresult without the 8(1 — zi)8(l — zi) 
part can be computed from our formula given in eqn.©. We can also obtain a result to N 4 LO 
order where we can predict partial sv contributions containing everything except the terms in 
<Do(zi)&(l-Zj),2>o{zi)2>o{zj),0i(zi)&{l-Zj) and 8(1 -zi)S(l -zi) for the coefficient A^ (4) . 
The results are identical to those given in the Appendix B of [36] for n\ = yep = M 2 . The convo- 
lutions of distributions of the form £>i(zj) <8> £> m (zj) for any arbitrary /,m can be done using the 
general formulae given in [30] so we obtain A s Jj® for i = 1,...,4. The differential cross sections 




Figure 1: Rapidity distributions for Z boson production at the LHC, and their /j = /jr (left panel) 
and /j = fj F (right panel) scale dependence (with M|/2 < /j 2 < 2M|). The abbreviation "pSV" 
means partial-soft-plus-virtual. 

can be expanded in powers of the strong coupling constant as 

d 2 ^ " . d 2 c J >^ 

dq 2 dy ~q s dq 2 dy 

We split the partonic cross section into hard and sv parts: 

dV'W _ d 2 G hard ' y 'W d 2 G sv > J >® 

dq 2 dy dq 2 dy dq 2 dy ' 
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Figure 2: Rapidity distributions for W + boson production at the LHC, and their h = hr (left panel) 
and /j = /jf (right panel) scale dependence (with M 2 v /2 < /j 2 < 2M(y). The abbreviation "pSV" 
means partial-soft-plus-virtual. 



2S — dq 2 dy — ~~ ^ ^ SM {^fi^^^A) + D qg '^(A^A^f) 

+D S 8 ^\AAti)) ■ (30) 
The SM coefficients D S a f' {i) (x° v x° 2 ^f) 

can be found in [6,46]. The sv part of the partonic cross 

section can be expressed as 

j2 a sv ,/,(() fl ri 

2S dq 2 dy = £- gSM Jo dXl Jo dx ^Kb^ x ^F) 

x £ dzi f Q dz 2 5(x? - xm ) 8(x§ - X2zi) A s jf (zi , z 2 , ^ 2 ,A/|,A/i) • (31) 

The coefficients A^ v a ^(zi,z2,<? 2 ,^,^) are presented in the Appendix B of [36], with the normal- 
isation A™J?\zi,Z2,q 2 ,Li F ,vl) = 5(1 -zi)5(l -zi). The functions g/ M are given by 

Mz _ 4oc 2 Q 2 , 4a< 7 2 r z ^ +r / v 2 , 2 x 

^ 3^ + Mz ((, ? 2_ M 2 ) 2 +M 2 r 2 )c 2 5 2^^ + W J 
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2a 2 (1 -4*2) (42 _ M 2 z) 



Qqg, 



'q > 



3(( ? 2-Ml)2+M|rl)c24 



(32) 



M w _ 



Oq r w ^iv, 




(33) 



((■ 



q 2 -M 2 v ) 2 +M 2 v T 2 w ^M W S 2 



We now give results for Z and W ± production by choosing q 2 = M| and q 2 = M w respectively. At 
these points the functions above contain the standard electro-weak constants which can be found 
in [6, 44] and the standard CKM matrix elements Vjj. We present our results as differential cross 
sections in rapidity for these fixed q 2 values and to compare with other authors plot d 2 o/dqdy 
where q = Mz or q = Mw respectively. 

We choose the center-of-mass energy to be \/S =14 TeV for the LHC and \/S =1.96 TeV for 
the Tevatron. The Z boson mass is taken to be Mz = 91.19 GeV and the width is 2.50 GeV. The 
corresponding values for the W boson are Mw = 80.43 and 2.12 GeV respectively. The strong cou- 
pling constant 0C v (/i|) is evolved using the 4-loop RG equations depending on the order in which the 
cross section is evaluated. We choose c£°(Af z ) = 0.130, oc^ LO (M z ) = 0.1 19, af LO (M z ) = 0.115 
and af LO (M z ) = 0.113 for i > 2. The set MRST 2001 LO is used for leading order, MRST2001 
NLO for NLO and MRST 2002 NNLO for N'LO with i > 1 [55, 56]. By choosing these parton 
densities we can compare our results with those of other authors (see later). We use a = 1/128 
for the electromagnetic fine structure constant, sin 2 0iy = 0.2314 for the weak mixing angle and 
cos 2 Be = 0.975 for the Cabibbo angle. In fig. 1 we plot the rapidity distributions for the Z-boson 
at the LHC in LO (dotted lines), NLO (solid lines), N 2 LO(SV only, dot-dashed lines) and N 3 LO 
(pSV only, small-dashed lines). Note that we have not plotted the partial sv N 4 LO contributions 
and we have only plotted the curves above 20 pb/GeV to magnify the central rapidity region. There 
are two panels in this plot and two curves in each panel since we show the scale variations by vary- 
ing the mass factorization scale fi F in the parton densities and the mass renormalization scale ijlr 
in the coefficient functions. Therefore we plot the curves at fixed /up = M| but with nj, = Af|/2 
and n 2 R = 2M| in the left panel and fixed fi\ = M§ but with f? F = M§/2 and i? F = 2M| in the right 
panel. Our results in the left panel show that there is only a tiny dependence on the jur for fixed 
/uf (here the LO result has no variation). In the right panel we see that /Uf dependence for fixed /jr 
decreases as we go from LO to NLO, N 2 LO and N 3 LO respectively. The N 3 LO band lies within 
the N 2 LO band and both are within the bands for the NLO results. We also notice that the lower 
curves for the N 2 LO and N 3 LO results fall on top of each other. The actual numbers are different 
but so close that one cannot see this from the plot. These results demonstrate that the perturbation 
series for the rapidity distribution converges very nicely at the LHC energy. 

In fig. 2 we plot the rapidity distributions for the W + -boson at the LHC in LO (dotted lines), 
NLO (solid lines), N 2 LO(SV only, dot-dashed lines) and N 3 LO (pSV only, small-dashed lines). 
We have only plotted the curves above 200 pb/GeV to magnify the central rapidity region. Again 
there are two panels in this plot and two curves in each panel since we show the scale variations by 
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Figure 3: Rapidity distributions for Z boson production at the Tevatron, and their /j = /jr (left 
panel) and n = hf (right panel) scale dependence (with A/f/2 < i? < 2M|). The abbreviation 
"pSV" means partial-soft-plus-virtual. 

varying the mass factorization scale /jf in the parton densities and the mass renormalization scale 
fXR in the coefficient functons. Therefore we plot the curves at fixed /up = but with fi R = M^/2 
and jd 2 R = 2M\ / in the left panel and fixed /j r = but with /j f = /2 and fjp = 2M^ in the right 
panel. Our results in the left panel show that there is only a tiny dependence on the /jr for fixed 
Hf (here the LO result has no variation). In the right panel we see that /jf dependence for fixed 
Hr decreases as we go from LO to NLO, N 2 LO and N 3 LO respectively. However the N 3 LO band 
is slightly below the N 2 LO band near y = even though both are within the bands for the NLO 
results, which is probably caused by the fact that we only have a partial soft-plus-virtual N 3 LO 
result. We notice again that the lower curves for the N 2 LO and N 3 LO results in the right panel fall 
on top of each other. The actual numbers are different but so close that one cannot see this from 
the plot. Both plots indicate that the perturbation series is rapidly converging at the LHC energy. 

In figures 3 and 4 we repeat these plots for the Tevatron energy and the same scale choices as 
above. The left panel in Fig. 3 shows excellent convergence of our results. Note that we only plot 
our results above 5 pb/GeV to magnify the central rapidity region. However the right panel shows 
that the bands for the N 3 LO result are wider than those for the N 2 LO result and both are wider than 
those for the NLO result. This can have two reasons. One is that our results are only sv or partial 
sv. The other is that since the Tevatron is an antiproton-proton collider different combinations of 
parton densities are involved as we increase the order of the perturbation series. We discuss this 
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Figure 4: Rapidity distributions for W + boson production at the Tevatron, and their fi = /jr (left 
panel) and /j = fif (right panel) scale dependence (with Mjy/2 < fj 2 < 2M^). The abbreviation 
"pSV" means partial-soft-plus-virtual. 

again below. We notice again that the lower curves for the N 2 LO and N 3 LO results in the right 
panel fall on top of each other. The actual numbers are different but so close that one cannot see 
this from the plot. 

Figure 4 shows results for W + production above 40 pb/GeV to concentrate on the central 
rapidity region. The curves in the left panel show excellent convergence of the perturbation series. 
The right panel again shows that the bands for the N 3 LO result are wider than those for the N 2 LO 
result and both are wider than those for the NLO result. Again we believe that this is caused by a 
combination of the reasons above and comment on it below. We notice again that the lower curves 
for the N 2 LO and N 3 LO results in the right panel fall on top of each other. The actual numbers are 
different but so close that one cannot see this from the plot. 

Note that the asymmetry about y = in the rapidity plots for W + production at the Tevatron in 
Fig.4 has basically disappeared at the LHC energy (see Fig. 2). 

We have checked our results in two ways. First by comparing our curves with similar plots 
for the rapidity distributions in [12]. Their computer program for the rapidity distributions has 
the exact LO result, the exact NLO result and a sv approximation for the N 2 LO result. We agree 
with their numbers when we choose their values for the electroweak parameters and their parton 
densities. Second we have also checked our results against those in [32], where the exact N 2 LO 
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rapidity distributions for the Z and W bosons are calculated. Their paper contains plots for the 
exact LO, the exact NLO and the exact N 2 LO results in pQCD. We have run their computer code to 
compare their results against ours. Our sv approximation agrees very well with their N 2 LO results 
in the case when ftp = /jr = Mj/2 for the Z and W ± bosons. However we get a slightly wider 
band when we vary the scales in the N 2 LO case, since we only have a sv approximation. From this 
comparison we can see that the wider bands we observe in Figures 3 and 4 are basically due to the 
sv and partial sv nature of our higher order results. We can probably reduce the width of the bands 
in the N 3 LO case by calculating the missing pieces of the partial sv result. However the central 
value of our partial sv N 3 LO result is very small when compared with the N 2 LO result indicating 
that the perturbation series continues to converge rather rapidly. Hence this is good news for the 
LHC experimenters who plan to calibrate the Atlas and CMS detectors by measuring the rapidity 
and transverse momentum distributions of Z and W ± bosons. 

To summarise, we have systematically studied higher order sv corrections to rapidity differ- 
ential distributions for Z and W ± boson production. We have used Sudakov resummation of soft 
gluons to calculate these processes. The resummation of soft gluons has been achieved using renor- 
malisation group invariance and the factorisation property of the observable that is considered here 
(the rapidity). Using the available information on the form factors, the DGLAP kernels and lower 
order results we have obtained compact expressions for the resummation of soft gluons for the 
rapidity distributions of Z and W ± bosons. Using these we have computed sv rapidity distributions 
exactly at N 2 LO and partially at N 3 LO. We have presented the numerical impact of these results. 
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